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We determine the phase diagram and chiral condensate for lattice QCD with two flavors of twisted- 
mass fermions in the presence of nondegenerate up and down quarks, discretization errors and a 
nonzero value of Oqcd- Although such a theory has a complex action and cannot, at present, be 
simulated, the results are needed to understand how to tune to maximal twist in the presence of 
electromagnetism, a topic discussed in a companion paper. We find that, in general, the only phase 
structure is a first-order transition of finite length. Pion masses are nonvanishing throughout the 
phase plane except at the endpoints of the first-order line. Only for extremal values of the twist 
angle and Oqcd (cj = 0 or 7r/2 and Oqcd = 0 or n) are there second-order transitions. 


I. INTRODUCTION 


In Refs. [1- and [2] we determined the phase structure 
of two-flavor lattice QCD with Wilson and twisted-mass 
fermions at nonvanishing lattice spacing in the presence 
of the two physical sources of isospin breaking: nondegen¬ 
erate up and down quarks and electromagnetism. These 
results are relevant for present simulations because 0(a 2 ) 
discretization effects are of comparable size to those from 
isospin breaking (Here a is the lattice spacing). Dis¬ 
cretization effects can thus significantly distort the chiral 
condensate and lead to unphysical phase transitions. 

In particular, the CP-violating phase found by Dashen 
in the continuum [3] can be enlarged by discretization 
effects, and for large enough a can include the point with 
physical quark masses. 

In Ref. [2 we found that the inclusion of electromag¬ 
netism along with twisting requires one to consider, at an 
intermediate stage, a lattice theory that has, in addition 
to isospin breaking, a nonvanishing value of Oqcd- We 
repeat the explanation of this result below. The purpose 
of the present note is to study the properties of this ex¬ 
tended theory, providing results that are used in Ref. |2] 
to tune to the physical value, Oqcd = 0- We also present 
some results not needed in Ref. |2] so as to provide a com¬ 
plete picture of the parameter dependence of the phase 
structure. 

Our analysis is carried out using SU(2) chiral per¬ 
turbation theory (XPT). Previous work has considered 
this theory at nonvanishing Oqcd in the continuum. In 
particular, Refs. [3] and [5] have studied the theory at 
Oqcd = tt, elucidating the Dashen phase. In addition, 
Refs. [6] and [7] considered this and related theories for 
arbitrary Oqcd in the small-volume regime, where one 
can use the methods of random matrix theory. In all 
these theories they find a Dashen phase. What our study 
adds to this previous work is the inclusion of the effects 
of working at nonvanishing lattice spacing. 


* e-mail: dhorkel@uw.edu 
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We begin by recalling the essential features of 
SU(2) XPT including discretization effects, nondegener¬ 
ate quarks, and electromagnetism, in the power count¬ 
ing we use in Refs, [l] and [2]. In this power count¬ 
ing, effects proportional the average light quark mass, 
m q = (m w + m^)/2, are assumed comparable to those 
quadratic in lattice spacing^ and to those proportional 
to c^em? *-e., m q ~ cl 2 ~ Q!em- We also include in the 
leading-order Lagrangian the dominant term introduced 
by nondegeneracy, which is proportional to e 2 , where 
e q = (m u — my)/ 2. We work to leading order in this 
combined power counting, so that loop effects need not 
be considered. Ignoring electromagnetism for now, the 
Lagrangian is then 

C x = f Atr[d^d^ t]+V x (1) 

V x = -L tr(x t S + ER) - Wqtr(A t X + eU)] 2 

+ L [tr(x t S _ S t x)] 2. (2) 

Here E G SU(2) is the chiral field, / « 92 MeV and 
are the continuum leading order low energy coefficients 
(LECs), and A = 2Woal is a spurion field, with Wo and 
W' LECs introduced by discretization errors. The quark 
mass matrix, M, is contained in the convenient quantity 
X = 2 HqMq Matching physical quantities in continuum 
SU(2) and SU(3) XPT one finds 


/ 2 

8B 0 m s 


(3) 


where m s is the strange quark mass 0. Thus we know 
that f'j is positive. 


1 Terms linear in a , if present, can, in the pion sector, be absorbed 
into an additive shift in the quark mass, so that the leading 
discretization effects relevant for the phase structure are propor¬ 
tional to a 2 [8]. 

2 The detailed relationship of the masses in M to the bare lattice 
quark masses is explained in Appendix A of Ref. [2]. We also 
note that M contains only the LR projection of the full mass 
matrix. 
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In the continuum, the leading term induced in the chi¬ 
ral potential by electromagnetism is that due to one- 
photon exchange between electromagnetic currents m 

m 


Vem = -^c E m tr(Er 3 i: t r 3 ), (4) 

where cem is proportional to «em and is known to be 
positive [12.. Electromagnetism also contributes to mass 
renormalization, but this is implicitly included by our use 
of renormalized masses in the quark mass matrix M. Al¬ 
though the quark masses depend on the renormalization 
scheme chosen, this dependence is canceled by that of 
the prefactor Bq, so that the product x is independent 
of renormalization scheme and scale. 

On the lattice, with Wilson or twisted-mass fermions, 
the inclusion of electromagnetism leads to additional is¬ 
sues. The first of these concerns the direction of the 
twist. Quark nondegeneracy picks out the 73 direction in 
isospin space. In the absence of electromagnetism, one 
can twist in an orthogonal direction, e.g. ti, and this 
choice leads to a real lattice fermion determinant m ■ 
However, such a twist leads leads to an electromagnetic 
current that includes an axial component when written 
in terms of bare quarks. This current cannot be coupled 
in a gauge-invariant way to the electromagnetic field in 
a lattice theory since it is not conserved [ 2 ]. 

To include both electromagnetism and nondegeneracy 
on the lattice, one is thus forced to twist in the 73 direc¬ 
tion mug. This, however, leads to a complex lattice 
fermion determinant m, making the theory challenging 
to simulate^] An intuitive way of understanding why the 
action is complex is to note that, in the continuum, when 
the twist angle is cj, the fermion mass term is given by 

TpirriqCu + T 3 e q c u + ii5T 3 m q s u + i^e q s u )ij} , (5) 


where is an isodoublet, c u = coscj, and = sincj. 
By construction, in the continuum a nonsinglet axial ro¬ 
tation (z.e., a twist) can return the mass matrix to its 
standard form m q -\-r^e q . However, on the lattice, such a 
rotation is not a symmetry. Crudely speaking, the lattice 
theory with mass matrix <© corresponds to a continuum 
theory in which the coefficients of the four terms are dif¬ 
ferently renormalized. In such a theory the mass terms 
involving 75 cannot both be rotated away, and thus the 
theory has a nonzero ©qcd- As is well known, this leads 
to a complex fermion determinant. The only redeem¬ 
ing feature is that, if one could tune the lattice quark 
mass matrix such that it took the form of Eq. (|5| in the 
continuum limit, then the imaginary part of the fermion 
determinant would vanish in this limit. 


3 This is avoided in Refs. m and m by working to linear order 
in perturbation theory about the isospin symmetric theory. 


The situation is not this simple, however, because 
of the second issue induced by the inclusion of elec¬ 
tromagnetism in the lattice theory. This is the pres¬ 
ence of independent additive renormalizations of the up 
and down bare untwisted quark masses proportional to 
c^em/cl. Since in our power counting m ~ a 2 ~ cuem , 
these renormalizations dominate over the leading order 
terms described above and collected in Eqs. © and ©>• 
They must be tuned away by applying nonperturbative 
conditions to determine, independently, the two critical 
masses]^] One of the results of Ref. [ 2 ] was a demonstra¬ 
tion that the tuning method used in Ref. m does not 
work in general^] The method provided only a single con¬ 
dition, while two are needed. The key point for present 
purposes is that, with the untwisted parts of the quark 
masses “detuned”, the theory one is studying has, even 
in the continuum limit, a nonvanishing value of ©qcd- 
Thus, to come up with a second condition that will set 
©qcd = 0 (in the continuum limit) one must understand 
the properties of the detuned theory. This is the purpose 
of the present analysis. 

To understand why detuning leads to nonzero ©qcd, it 
is instructive to write out the renormalized mass matrix 
M in a detuned, twisted theory. It is convenient to work 
with y rather than M, since this is what enters the chiral 
Lagrangian. The form is 


X = 


(fh + */i u 0 

V 0 - ifid) 


(6) 


where the “hat” on a mass indicates multiplication by 
2E>o- In particular, = 2£>ora^, with the renor¬ 
malized untwisted or “Wilson” part of the up-quark 
mass, while Jl u = 2£?o/i u , with /jl u the twisted part of 
the up-quark mass. Similar notation holds for the down- 
quark masses. The superscript W distinguishes the un¬ 
twisted masses from the full physical masses, which are 
given, for example, by m 2 = (mj ^) 2 + /i 2 . The two 
twisted masses in (J 6 | have opposite overall signs because 
twisting involves 7-3. In this notation, tuning the un¬ 
twisted parts of both masses to their critical values means 
tuning both and fhY to zero. As long as fi u and fid 
have opposite signs this corresponds to tuning to maxi¬ 
mal twist. 

We can rewrite the mass matrix of Eq. (| 6 | in terms of 
the average physical quark mass m q and the nondegen- 


4 There is, in addition, the standard additive renormalization pro¬ 
portional to 1/a (times powers of a s ) that is common to both 
quarks. The nonperturbative conditions that remove the Q!em/® 
shifts will also remove the larger 1/a shifts. The point here is 
that the smaller (but still divergent) electromagnetic renormal¬ 
izations imply the need for two conditions, rather than one. 

5 The method, based on introducing unphysical valence quarks, 
works only in the electroquenched approximation, in which sea 
quarks are kept neutral and degenerate. It fails once the sea 
quarks are charged. The second method proposed in Ref. mi. 
and the method used in Ref. E2, do not suffer from the same 
problem, because they tune using physical quantities. 
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eracy e q : 


X = 


(fhg + £q)e^ v+w ) 


(fhg - e g )e i( ^ w) 


(7) 


Here m q = 2 B 0 m q , e q = 2 B 0 e q , 

tan(f + uS) = and tan(< p — iS) = —-^^. ( 8 ) 

m u m d 

We observe from Eq. 0 that c o is the twist angle, while 
p, being the overall phase of the mass matrix, is propor¬ 
tional to 0QCD- 


= 


Qqcd 

2 


(9) 


Thus having a general, detuned mass matrix corresponds 
to working at nonzero Oqcd- Tuning to the critical val¬ 
ues of the untwisted quark masses corresponds to setting 
uj = 7r/2 and p = 0, i.e., tuning to maximal twist with 
vanishing 0qcd- 

In summary, the dominant effect of including elec¬ 
tromagnetism in a theory with Wilson or twisted-mass 
fermions is mass renormalization. For nonvanishing 
twist, this implies that one must work at nonvanishing 
0qcd in order to tune to 0qcd = 0. While this will 
be challenging for simulations, it is straightforward to 
study this theory in XPT. Working at leading order in 
our power counting, one has simply to find the minima 
of the potential that is composed of the terms given in 
Eqs. ([ 2 ]) and Q. 


When determining the expectation value of the chiral 
field, it is convenient to parametrize it relative to the 
twist it would obtain were <p = a = 0 : 


<E> 


_ e iujr 3 /2^iOh-T e icdr 3 /2 


( 10 ) 


The full potential V = V x + Vem then becomes 


V 


f 2 


fh q cos 9 cos cp + nse q sin 0 sin cp 

+ C£ ( ns€ q sin 6 cos cp — fh q cos 6 sin p ) 2 
+ w' (ns sin 9 sin uj — cos 9 cos uj ) 2 
+ c E M {n\ + (1 - n\) cos# 2 ) , 


( 11 ) 


up to an irrelevant overall constant. Here we have intro¬ 
duced 


An immediate conclusion from this analysis is that the 
remaining explicit effect of electromagnetism, namely the 
cem term, is simply a constant for ns = 1. It therefore 
does not effect the minimization of the potential, and 
thus has no impact on the phase structure. Physically 
this is because the condensate lies in the neutral pion 
direction. The only effect of this term is togive an overall 
positive shift in the charged pion masses jj 


The remainder of this paper is organized as follows. In 
Sec. [IT| we determine the phase diagram in the fh q — e q 
plane. We do so in stages, beginning by elucidating the 
symmetries of the potential ( 11 ), then working out the 
phase diagram in the continuum, next adding in dis¬ 
cretization effects for the extremal cases where uo = 0 
and 7r/2, and finally considering the most general choices 
of parameters. We then return, in Sec. m to the original 
motivation for the present work, namely the determina¬ 
tion of a condition such that, in the presence of electro¬ 
magnetism, maximal twist at 0qcd = 0 can be achieved 


in a physical phase. We conclude in Sec. IV 


II. DETERMINATION OF PHASE DIAGRAM 


A. Symmetries of the phase diagram 


Before entering into detailed calculations we collect 
some general results that follow from the form of the 
potential, Eq. 

First we note that, without loss of generality, we need 
only consider uj and p in the range 

0 < uj, p < 7r/2 , (13) 

as long as we consider the full fh q — e q plane. This is 
because V is invariant under each of the following four 
transformations 


(i) {UJ —>• UJ + 7r} , (14) 

(ii) —y — uj , 6 —y — 0 , 6 q —y —, (15) 

(in) {if -> - f, e q -e q } , (16) 

(iv) {f —> ip + 7r, m q —> — m q , e q —> —e q } . (17) 


In the following, we refer to the endpoints of the range 
(13) as the “extremal” values of uj and p, while values 
within the range are called “nonextremal”. 

In addition, V is invariant under 



MW'Wfia 2 

T 2 


and w' = 


( 12 ) 


Given that £7 and cem are both positive, the potential 
is always minimized with the condensate aligned in the 
t 3 direction, i.e. n = (0,0, ± 1 ). Without loss of gener¬ 
ality we can set h = ( 0 , 0 , 1 ) and absorb any sign into 0 . 
The main task in the following is the determination of 
the values of 6 which minimize V as the parameters are 
varied. 


(v) {uj —y — —uj, p ^y ——p, m q O e q , 0 — — 0 } . 

(18) 

This implies relations between the phase transition lines 
for different values of the parameters. For example, the 


In light of these considerations,we drop the cem term in the 
subsequent discussion of minimization of the potential. 
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phase diagram for uj = p = 0 is related to that for uj = 
p = 7r/2 by a reflection in the diagonal line fh q = e q . 

The final invariance that plays a role in the following 
is 

7T 

(vi) \uj —y -cj, £q —y — c q , w —y —w , # —y —#j-. 

(19) 

This relates the phase diagrams with opposite signs of 


B. Continuum XPT with nonzero ©qcd 


The absence of a transition for nonextremal p contin¬ 
ues to hold when the eg term is restored. This can be 
understood as due to the lack of a Z 2 symmetry in the 
potential. It is the presence of a Z 2 symmetry for ex¬ 
tremal p (under which # —>> # + it) that, when broken 
by the vacuum, leads to a second-order transition. The 
upshot is that the extremal phase diagrams of Fig. [l] are 
replaced by blank diagrams with no transitions, aside 
from the singular point at the origin. 

To show a concrete example of this, we consider 
p = 7r/4. Using the parametrization fh q = rcosa, 
e q = rsina and r 2 = ra 2 + e 2 . The potential is then 


In this section we examine the phase structure in the 
continuum. Without discretization effects, the twist an¬ 
gle is redundant and has no effect on the phase diagram. 
This is manifest in the basis used in Eq. ( pTj ), where with 
w' oc a 2 = 0 there is no dependence on cj. 

We begin by recalling results for the extremal cases 
p = 0 and 7r/2, corresponding to Oqcd = 0 and 7 r. These 
have real and imaginary quark masses, respectively. 

The physical case, p = 0, has been described exten¬ 
sively in the literature mmm- The phase diagram 
is shown in Fig. [lja). There is a second-order transi¬ 
tion between the standard continuum phase and the CP- 
violating Dashen phase, lying along fh q = =L2qc 2 . In the 
Dashen phase, the potential has two degenerate minima, 
both having 


cos # = 


2 cge q ’ 


( 20 ) 


and differing in the sign of #. 

The case of ©qcd = n (p = 7t/2) was first described 
by Smilga [4]. The potential has the same form as for 
p = 0, except that fh q and cos# are exchanged with 
e q and sin#, respectively. This implies that the phase 
diagram has the same form as for p = 0, except that it 
is reflected in the fh q = e q line, as shown in Fig 


•0b). 

(since 


This is an example of the symmetry ( [18] ) at work 
the change in uj is irrelevant in the continuum). There 
is thus a second-order transition to a Dashen-like phase 
along the lines e q = ±2 cgfh? q . There is again a two-fold 
degeneracy within this phase. 


For nonextremal p (0 < p < 7t/2) the potential (11) 
cannot be minimized exactly and it is instructive to look 
at some simple limits. 

First we drop the 0(m 2 ) eg term. In the extremal cases 
this means that the width of the shaded (pink) phases 
shrinks to zero, so that there is a first-order transition 
along the entire fh q = 0 line for p = 0 and along the 
e q = 0 line for p = 7t/2. By contrast, for nonextremal p, 
there are no transitions. The potential is minimized at, 




a) + ft cos 2 (# + a)] , 


( 22 ) 


where ft = qr/ y/2 can be treated as small in our power 
counting. The minima occur when 

0 = sin (# — a) + 2ft sin [2(# + a)] . (23) 

Expanding in powers of ft about the leading-order solu¬ 
tion, # = we find 

6 = a — 2ft sin (4a) + 0(ft 2 ). (24) 


The presence of only a single solution indicates the ab¬ 
sence of a Dashen-like phase. We have investigated this 
numerically for other values of p and found that there 
are no phase transitions for any nonextremal p. 

To see how the degeneracy of the Dashen-like phase is 
broken for nonextremal p, consider the potential along 
the fhq = 0 axis: 



e q sin # sin p + Qe 2 sin 2 # cos 2 p . 


(25) 


For p = 0, one finds (since eg > 0) that there are degen¬ 
erate minima at sin# = =bl. This corresponds to moving 
from the origin in Fig. flVa) along the e q axis and thus 
lying in the (shaded pink) Dashen phase. Turning on a 
nonzero p, the potential is still extremized at | sin # | = 1, 
but the two extrema are no longer degenerate 


V(sin# = ±1) 

T 2 


sin p + Qe 2 cos 2 p . 


(26) 


Thus there is a unique minimum, such that sin # = 1 for 
e q > 0 and sin# = — 1 for e q < 0 (assuming a positive p). 
There thus can be no Dashen-like phase. 


C. Discretization effects at nonzero ©qcd for 
extremal uj 


tan # = tan p , (21) 

m q 

and changes continuously as one moves through the phase 
diagram, except when passing through the origin. 


We now turn on discretization errors by considering 
non vanishing w'. Just as in the continuum, the phase 
diagram is easiest to determine for extremal p. The case 
of uo = p = 0 (untwisted fermions with ©qcd = 0) has 
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rh q rh q 




FIG. 1: Continuum phase diagram for (a) p = ©qcd/2 = 0 and (b) p = n/ 2 . Shaded (pink) regions have varying 
values of the vacuum angle #, as indicated in the figures. Unshaded regions have constant #. The neutral pion mass 

vanishes along the phase transition lines. 


long been studied, and it has been shown that there are 
two distinct scenarios depending on the sign of w'\ the 
so-called Aoki scenario for w' < 0, and the first-order sce¬ 
nario for w' > 0 mm- The resulting phase diagrams 
are shown in Fig. [2j and should be compared to the con¬ 
tinuum diagram of Fig. [lja). For w' < 0, the Dashen 
phase, in which # is degenerate, expands vertically so 
as to include the origin. The CP violating phase along 
the fh q axis is typically called the Aoki phase, so we call 
the extended CP violating region the Aoki-Dashen phase. 
This situation is shown in Fig.[2ja). For w' > 0, the ver¬ 
tical width of the continuum Dashen phase is reduced, 
and there is a segment of first-order transition along the 
e q axis, as shown in Fig. [2|b). In both scenarios, within 
the Aoki-Dashen phases the potential is minimized by 


cos 0 = 


frig 

2 (c e e q - w’) ’ 


(27) 


so that again there are two degenerate vacua with oppo¬ 
site signs of #. 

We next consider p = n /2 while holding uj = 0 (i.e. 
Wilson fermions at ©qcd = tt). This has not been pre¬ 
viously studied in the presence of lattice artifacts. As 
described above, for the continuum terms in V, chang¬ 
ing p from 0 to 7r/2 has the effect of interchanging fh q 
and cos# with e q and sin#, respectively. Since the w' 
term can be rewritten as w' cos 2 # = w'( 1 — sin 2 #), the 
same interchanges hold for w' ^ 0 as long as one flips 
the sign of w f . Up to some unimportant sign flips, this 
is an example of the general transformation obtained by 


combining Eqs. (1181) and (19): 


(' vii) {p^ --p, fa q 


w' — —w\ #—)►#——}. (28) 


The implication is that the phase diagrams for p = 7r/2 
are obtained from those of Fig. [2] by rotating 90° coun¬ 
terclockwise, and interchanging the w' < 0 and w' > 0 
scenarios. The positions of the resulting transitions are 
shown schematically in Fig. [3j 

The twist angle is no longer redundant when w' ^ 0, 
entering the w' term in Eq. (ITlJ) as w' cos 2 (# + uS). Thus 
changing uj from 0 to 7r/2 has the effect of flipping the 
sign of w': w' cos 2 (# + ir/2) = w' sin 2 # = w'( 1 — cos 2 #). 
This is an example of the general transformation ( p~9] ) . It 
implies that the phase diagrams for maximal twist can be 
simply obtained from those without twist. The situation 
is summarized in Fig. [3] 

In the remainder of this subsection we keep uj at one of 
the extremal values but allow p to take on nonextremal 
values. We recall that in the continuum, the phase di¬ 
agram with such parameters has no phase transitions. 
This turns out not to be the case when w' ^ 0. Exam¬ 
ples of the results we find are shown in Fig. [4j 

We begin with uo = 0 and nonextremal p, and work 
in the w' > 0 scenario. We find that there is a first- 
order transition along a finite segment of the e q axis, 
across which # changes discontinuously. The length of the 
segment depends on p. As p approaches zero jin which 
limit one obtains the phase diagram of Fig. [2|b)] the 
first-order segment asymptotes to precisely the first order 
transition line shown in Fig.[2|b), with end points qc 2 = 
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rh q m q 




FIG. 2: Phase diagrams from Ref. PQ including effects of discretization for uj = p = 0: (a) Aoki scenario (w' < 0) 
and (b) first-order scenario (w' >0).. The expression for 0 in the shaded (pink) region in (a) also holds in (b). The 
boundary of the shaded regions are second-order transition lines, along which the neutral pion mass vanishes. The 
(yellow) solid line running along the e q axis between the shaded regions is a first-order transition. 



w'<0 w=0 <p=0 and w'>0 w-n!2 <p=0 
w‘>0 (x>- 0 <p=0 and w'<0 (jj-n!2 <p=0 
w‘>0 w=0 <p=nl2 and w'<0 a)=nI2 <p=nI2 
w'<0 w=0 <p=nI2 and w >0 w=n!2 (p=n!2 


FIG. 3: Schematic positions of phase boundaries for extremal choices of uj and cp for both w' < 0 and w' > 0 
scenarios. Dashed lines indicate second order transitions, solid lines indicate first-order transitions. Results for two 

of these parameter choices also appear in Fig. [2] 


Trek We stress again that, for non vanishing (/?, there 
are no regions of Aoki-Dashen phase. As (p increases, the 
first-order segment reduces in length, until, as ip —>> 7r/2, 
it approaches the width of the Aoki phase that appears 
at p = 7r/2, i.e. with end points e q = ±2 w'. [Recall 
that the phase diagram at ip = tt /2 is given by Fig. |2][a) 
rotated by 90°; see also Fig. [3j] The first-order segment 
at the halfway point, p = 7r/4, is shown as the horizontal 
solid (red) line in Fig. [4j 


The length of the segment can be obtained analytically 
for all p. To do so, one extremizes the potential after 
setting rh q = 0. The global minimum lies at 


sin# = 


— e q sin(^ 

2(q^ cos 2 p — w') 


(29) 


with the sign of cos 6 undetermined. As one passes 
through the transition line (by varying m g ), cos 0 changes 
sign, indicating a first-order transition. Solving for the 
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endpoints, where cos 0 = ±1, we find 

. „ . — sin p + \J sin 2 p + 16qk/ cos 2 p 

\^Q\ - A o 

4 C£ cos 2 p 


(30) 


This gives the results quoted above in the limits ip 
0,7r/2. 

The corresponding results for w' < 0 can be obtained 
from those just described for w' > 0 using the trans¬ 
formation of Eq. ( [28] ) . In words, to obtain the phase 
diagram for p = po and w' = w' 0 < 0, one takes the dia¬ 
gram with p = 7r/2 — po and w' = \w' 0 \ and rotates it by 
90° counterclockwise. This implies that the first-order 
transition line is now vertical. 

Similarly, one can obtain results for uj = 7r/2 from 
those at uj = 0 using the transformation of Eq. (18). 
Specifically to obtain the phase diagram for p = po at 
uj = 7r/2, one takes that with p = tt/2 — po and uj = 0 
and reflects it in the fh q = e q line. This implies that the 
first order line is vertical for w' > 0 and horizontal for 
w' < 0. An example of this result (for w' > 0) is shown 
by the vertical solid (purple) line in Fig. [4j 

Since there are no second-order phase transitions, 
the pion masses are nonvanishing throughout the phase 
plane, with the exception of the endpoint of the first- 
order transitions, where the mass of the neutral pion 
vanishes. 


D. Nonextremal uj and p 


Finally, we consider the most general choice of parame¬ 
ters: both uj and p nonextremal. Here, in most cases, we 
have to proceed numerically, although we can check the 
results analytically for the special case of uj = p = 7r/4. 

We have found in the previous subsection that, for ex¬ 
tremal uj but nonextremal p there is a first-order transi¬ 
tion line of finite length that is oriented either horizon¬ 
tally or vertically. For example, Fig. [4] shows that for 
p = 7r/4 and w' > 0, the transition line is horizontal 
at uj = 0 and vertical at uj = 7r/2. It is not surpris¬ 
ing, therefore, that for intermediate values of uj there is 
a first order transition line of finite length at an inter¬ 
mediate angle that interpolates between the horizontal 
and vertical limits. Examples for several intermediate 
values of uj for p = 7r/4, i r/8 and 3tt/S are shown in the 
figure. We observe that, aside from the special case of 
uj = p = 7r/4, the first order lines are “S-shaped” rather 
than straight. We also observe an example of an overlap¬ 
ping transition line (though of different lengths) for the 
parameter choices (uJ,p) = (7r/8,7r/4) and (7r/4,7r/8). 
We have not understood this overlapping analytically, 
and do not know if it is exact. Figure [4] also shows an ex¬ 
ample of the application of the symmetry of Eq. (18), 


which implies that the transition lines for parameters 
(uj,p) = (7r/8,7r/8) and (37r/8, 37 t/ 8) should be related 
by reflection in the diagonal fh q = e q line. 

For uj = p = 7r/4, we know from the symmetry of 
Eq. (18) that the transition line must be invariant under 


reflection in the fh q = e q line. Thus it must either lie 
along this line or be perpendicular to it, and in both 
cases it must be straight. It turns out that, for w' > 0 
it lies along the diagonal, as shown in Fig. [4] by the solid 
(blue) line. Given this information, it is straightforward 
to determine the end points analytically, and we find that 
they he at 


Tflq — € q — zb 


1 — \/l — 16C£W f 

4q 


(31) 


We note the curious result that these points lie at the 
junction of the boundaries of the Aoki-Dashen phases for 
extremal uj and p. 


III. MAXIMAL TWIST CONDITION 


The standard technique for tuning to maximal twist in 
the absence of isospin breaking is to enforce the vanishing 
of the “PCAC mass”. This determines the critical value 
of the untwisted component of the common quark mass. 
This critical value can then be used even in the presence 
of isospin breaking due to quark masses, i.e. when the 
twisted components of the up and down quark masses 
differ. This is no longer the case when electromagnetism 
is included, because, as explained in the Introduction, the 
critical masses for the up and down quarks differ. Setting 
the PC AC mass to zero is essentially a way of enforcing, 
in a particular correlation function, the restoration of 
SU(2) flavor and parity symmetries at nonzero lattice 
spacing. In the presence of electromagnetism, however, 
these symmetries are absent even in the continuum limit, 
so it makes no sense to enforce them. Thus one must use 
alternative methods to tune to maximal twist. 

In our companion paper, Ref. [2], we analyze a method 
for carrying out the tuning in the presence of electro¬ 
magnetism, proposed in Ref. HE This involves partial 
quenching, and our analysis is somewhat involved, but 
the details do not matter here. Our key finding is that 
the method fails to tune the untwisted components of 
the up and down quark masses to zero, as required for 
maximal twist, but rather only enforces a condition on 
the condensate: 

(S) = e^+^ T3 = e^/ 2 ^ 3 =► O'= 6 + u> = n /2 . (32) 


In the first equality we are simply using the definition of 
the phase angle #, given in Eq. (10), together with the 
result that n points in the 73 direction. This implies that 
the total t 3 rotation angle is 6 ' = 0 + cj, which is set by 
the condition of Ref. [15] to tt/2. 

The condition (32) is indeed consistent with the desired 
parameters, i.e. with uj = tt/2 and p = 0. To see this 
we note that, for these parameters, the phase diagrams 
are those of Fig. [ 2 ] except that the w' < 0 and w' > 0 
diagrams are interchanged Q Thus, as long as the physical 


7 This result is obtained by acting with the transformation ( |l9| 













— oj-O^-niA 

± t ^=7t/8 5 <p=jtM 
oj-nlA^-nlA 
¥ 0)=3nl8 I <p=nlA 

— w=7T/2 5 <p=7T/4 

* w=7r/8 5 <p=7r/8 
■ w=7t/4 5 ^=jt/8 

* w=3jt/8 5 <p=37t/8 


FIG. 4: Phase diagram including discretization effects (with w' > 0) for several values of uj and p. Solid lines are 
analytically determined first-order transition lines described in the text. Points represent the location of the 
numerically determined first-order transition lines. Dashed lines show the positions of the central second-order 
transition lines that arise at extremal values of uj and p as shown in Fig. [3| These are included to set the scale, since 
they depend on the values of w' and q. Results for w' < 0 can be obtained from these using the transformations of 

Eqs. (19) and (28). 


masses are such that one is in the unshaded region, i.e. as 
long as one avoids the Aoki-Dashen phases, one has that 
cos 0 = 1 and thus 6 = 0. This means that O' = uj = 7r/2, 
which satisfies Eq. (32). 

However, there is in general a one-dimensional family 
of solutions to Eq. (32), all having different values of uj 
and p. One can understand this intuitively as follows. 
Nonvanishing values of uj and of p both violate parity, 
and thus both lead to a nonzero twist of the conden¬ 
sate, i.e. a nonzero value of O '. For any choice of p, the 
desired value O' = 7r/2 can, in general, be obtained by a 
suitable value of uj. Thus there is a line in the uj — p plane 
along which the condition is satisfied. In order to tune 


on the uj = cp = 0 results that are actually shown in the figure. 
See also also Fig. [3] 


to the desired point on this line an additional condition 
is needed. 


It turns out to be easier to do the calculation using 
the parametrization of the mass matrix given in Eq. <§>• 
Here we fix the twisted components of the masses fi u and 
fid (ultimately to their physical values, namely 2F 0 mP hys 


and 2Hom^ hys , respectively) and vary the untwisted com¬ 
ponents ray and . This corresponds to what is done 
in actual simulations. The condition of Eq. (32) then 
forces the theory to lie on a line in the plane. 

Our aim is to determine this line and to find an addi¬ 
tional condition that picks out the desired point on the 


line, namely rajf = frfifi = 0. 
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In terms of the parametrization 0 the potential i 


is 


V_ 

'p 


rn 


w 


a' , fiu + fid • a , 
—— cos 0 H-sin 0 


+ Q 


l^u f^d /j/ Wl 

---cos 0 


W 


m 


w 


- sin O' 


+ w' cos 2 0 f . 


(33) 


In order for an extremum of this potential to lie at 0 ' = 
7r/2, it is simple to show that the untwisted masses must 
satisfy 


m 


w 


iW 


1 ~ cg(fi u - fi d ) 

1 + ctifiu - fi d ) 


= s. 


(34) 


i.e., the theory must lie along a straight line in the m— 
mY plane with slope s determined by the physical masses 
and eg. We can turn this into a constraint on cj and cp by 
equating the parametrizations of Eqs. 0 and ©. One 
finds fn^ d = fi u ,d cot (cj d= cp) , so that the allowed values 
of cu and cp satisfy 

s fi u cot (cj + (p) = fid cot (c o — (p). (35) 

Lines satisfying this equality are shown in Fig. [5] The 
desired point is at (p = 0, cj = 7 r/ 2 , but, as claimed above, 
solutions exist for all values of ip. 


U) 



FIG. 5: Values of cj and ip for which 0' = tt /2, using 
fi u « 2 Bom^ tlys and fid ~ 2 Bom d ciys , but with eg larger 
than the physical value so as to increase the curvature 
of the lines for the sake of clarity. 


The above considerations assume that the extrema at 
6 ' = 7 r /2 is a minimum of the potential. This can be 
determined by examining either the second derivative of 
the potential evaluated at the correct value of s and O’ = 
7r/2 or, equivalently, by checking that the neutral pion 


mass is nonnegative. The neutral pion mass along the 
line (l34|) is 


2 fiu~ffid 0 
m 7r° =—- 2c ^ 


l^u l^d 


+ 2 eg 


- w - w 

m u ~ m d 


-2 w'. 
(36) 


The sum of the first two terms is positive for physical 
parameters (since this is just the physical neutral pion 
mass-squared at this order in XPT, and higher order cor¬ 
rections are small). The third term is always positive and 
vanishes only at the point of maximal twist (as long as 
s 7 ^ 1, which is the case for physical parameters). The 
last term can be negative if re' > 0. Thus, if w' takes 
a large enough positive value, it can be that the point 
we are aiming to tune to does not lie at the minimum of 
the potential. This happens when the physical point lies 
inside the Aoki-Dashen phase. 

Assuming that this does not happen, we can ask 
what criterion can be used to tune to maximal twist 


along the lines satisfying Eq. (32). The criterion pro¬ 


posed in Ref. [ 2 | is simply to minimize the neutral pion 
mass, Eq. (36), since, as already noted, this occurs when 




= 0 . One can also minimize the charged pion 
mass, the expression for which is given in Ref. [ 2 ]. 


We close this section by making a connection with our 
results for the phase diagram for general cj and cp, ob¬ 
tained in Sec. [TT] In particular, we imagine that we have 
somehow tuned close to maximal twist, but that there 
is a small offset. Specifically wefixc<; = 7 r/ 2 -b( s ) and 
(p = aS with \S\ 1 and a ~ 0(1). This differs from 

(and is less realistic than) our analysis above where we 
fixed the twisted up and down masses. Nevertheless, this 
allows a valid theoretical exercise: with cj and (p fixed in 
this way, we determine the line in the fh q — e q plane that 
satisfies the tuning condition of Eq. ( [32] ) . 

Written in terms of the variables of Eq. 0 , the tuning 
condition becomes 

s(fh q + e q ) cos (cj + tp) = (fhq — e q ) cos(cj — cp) , (37) 

where we recall that s is given by Eq. ( [34] ) . For our fixed 
values of cj and cp, this equation can be converted into a 
result in the fh q — e q plane, using an expansion in powers 
of (5: 


a 




1 — 2 acgi q 


+ 0 ( 6 2 ). 


(38) 


Examples of this result (for both signs of w') are shown 
in Fig. [ 6 ] by the (green) solid lines. As noted above, this 
result is only valid if the pion mass-squared of Eq. (36) 
is positive or zero. Thus the line terminates at the point 
where m^o vanishes, which occurs only for w' > 0 . 

Also shown in the figures are the positions of the first- 
order lines, which have been determined numerically. We 
observe that, for w' < 0, the line along which Eq. (32) 


holds goes all the way to the origin, where it runs into the 
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first-order line. By contrast, for w' > 0 the endpoint of 
the tuned line is precisely the starting point of the first- 
order line. This is reasonable since it is the only position 
in the phase diagram where a pion is massless. In any 
case, we see that, even for nonextremal uo and ip, where 
there is only a first-order transition, the condition (32) 
cannot be maintained all the way to the origin. 


IV. CONCLUSIONS 

In this short note, we have determined the phase struc¬ 
ture of lattice QCD in the presence of isospin break¬ 
ing and a nonvanishing value of Oqcd- This is, for the 
present, a theoretical exercise, but one that was neces¬ 
sary in order to understand how to tune to maximal twist 
in the presence of electromagnetism, an analysis that was 
completed in our companion paper [2]. 

The results are also interesting in their own right. In 


particular, for generic (nonextremal) values of the twist 
angle and ©qcd, the continuum theory has no phase 
structure, while the lattice theory has a segment of first- 
order transition whose length is set by w' and is thus of 
0(a 1 2 3 4 5 6 7 8 9 ). 

We have kept in our analysis only the leading order 
terms arising from each type of symmetry breaking. A 
quantitative analysis would require the inclusion of all 
other terms of 0 {m 2 ) as well as those proportional to 
ma. Based on our work in Ref. [I], however, we do not 
expect these terms to lead to qualitative changes in the 
phase diagrams. 
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FIG. 6: Applying the tuning condition at fixed cj = 7t/2 + 0.1 and p = 0.05 for (a) w' < 0 and (b) w' > 0. The size 
of C£ and w' can be seen from the tick marks on the axes, which lie at fh q = ±2w' and e q = ±^\w'\ /q. The solid 
(green) line shows the result of applying the condition (32) as well as requiring that m^. 0 > 0. (Red) points show the 

locations of the numerically determined first-order transition lines. 










